X-ray Signatures of Circumnuclear Gas in AGN 
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Abstract. 

X-ray spectra of AGN are complex. X-ray absorption and emission features trace gas covering a 
wide range of column densities and ionization states. High resolution spectra show the absorbing gas 
to be outflowing, perhaps in the form of an accretion disk wind. The absorbing complex shapes the 
form of the X-ray spectrum while X-ray reverberation and absorption changes explain the spectral 
and timing behaviour of AGN. We discuss recent progress, highlighting some new results and 
reviewing the implications that can be drawn from the data. 
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INTRODUCTION 

Spectroscopic data from the X-ray, UV and optical regimes have established the pres- 
ence of significant columns of material outflowing from supermassive black holes in 
active galactic nuclei (AGN). Studying these outflows offers insight into the observed 
relationship between the black hole and its host galaxy [e.g. 1]. Given the short duty- 
cycle of accretion that seems to exist for AGN [2], epochs of significant nuclear outflow 
must occur for AGN to build the supermassive black holes observed in the Universe. 
The most promising models for X-ray outflows consider mass loss via an accretion disk 
wind [e.g. 3]. The X-ray band carries information about the material flow very close 
to the black hole, and offers important constraints on fundamental parameters of the 
system, such as the wind launch radius and acceleration mechanism. 

The spatial resolution of current X-ray instruments does not allow us to image the in- 
nermost regions of AGN directly. However, recent long X-ray observations, particularly 
those afforded by Suzaku and XMM-Newton have allowed a significant step forward us- 
ing spectroscopic and timing analyses. Here, we review our current understanding of 
AGN absorber systems, highlighting key contributions from the Suzaku satellite. 



X-RAY ABSORPTION IN AGN 

X-ray grating spectroscopy has allowed more than an order of magnitude improvement 
in spectral resolution at 1 keV, and a factor ~ 4 improvement at 6 keV, compared with 
CCD instruments. Accurate measurement of the energies of soft-band absorption lines 
such as those from O, Ne, Mg, Si and S, first traced relatively low ionization gas, out- 
flowing at hundreds to thousands of kms^^ [4, 5]. Such so-called 'warm absorbers' are 
common in AGN. The best-studied, local AGN reveal multiple zones of absorbing gas. 
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FIGURE 1. Unfolded Chandra HEG spectral data of MCG-6-30-15 from 2004. The spectral model is 
shown as a solid line. See text and Miller et al. [10] for details. 

whose signatures can be separated with grating data. Application of photo-ionization 
models can yield the ionization-parameter for the gas, and an estimate of the product 
of density and radial distance from the source. In principle, the gas density may be es- 
timated using sensitive triplet line ratios, however, in practice it has proved difficult to 
obtain interesting constraints. Absorber variability can also be used to constrain the gas 
density via consideration of the recombination timescale, although attempts to date have 
been limited by scant variability measurements. 

Estimates of the gas location vary. Some new work considers reverberation analysis 
within the X-ray band, allowing a constraint on the location of reprocessing gas out of 
the line of sight. A few cases studied to date place the reprocessor tens to hundreds of 
gravitational radii from the nuclear source [6]. Other methods include measurement of 
the gas ionization parameter that, combined with velocity constraints and consideration 
of reasonable limits on the gas geometry, yield radial estimates that range from within 
the optical broad-line-region (BLR) to outside the putative molecular torus [5]. 

An important development has been the discovery of deep K-shell absorption lines 
from H-like and He-like species of Fe [e.g. 4, 7, 8]. Suzaku has been effective at detection 
of these lines, owing to its high sensitivity at 6 keV. Importantly, Suzaku data yielded 
detections of Fexxv and FeXXVl in MCG-6-30-15 (Fig. 1), with outflow velocity 
~ 1800km s^^ [9, 10]. In general, measured line depths indicate the presence of very 
high columns of gas, A^h ~ 10^-^ — 10^^ atoms cm^^ in the line-of-sight, with velocities 
typically much higher than those detected in the soft X-ray regime [e.g. 11, 12, 8]. A 
sample study has found velocities up to ~ 0.3 c for these high-ionization absorbers [13] 
with some systems transporting significant mechanical energy from the nucleus to the 
host galaxy [14]. 

Systems detected to date span six orders of magnitude in ionization parameter, and 
five orders of magnitude in integrated column density. In summary, we have detected 
absorbing gas over the entire range of parameter-space to which current X-ray instru- 
ments are sensitive. 



SPECTRAL VARIABILITY 



In light of the obvious importance of complex absorption in MCG-6-30-15, Miller et al. 
[10] compiled data from Chandra, XMM-Newton and Suzaku to conduct the most com- 
prehensive broad-band analysis (0.5-45 keV) of a single AGN to date. Spectral vari- 
ations were decomposed using the SVD-PCA method, showing that spectral variabil- 
ity in MCG-6-30-15 could be fully described with a simple variations of one partial- 
covering zone of the complex X-ray absorber. The variable-covering absorber suggested 
for MCG-6-30-15 is of general interest, as it accounts for the spectral variations seen in 
a number of AGN [12, 10]. 

In the case of a few sources, such as NGC 3516 [8], the data require an additional com- 
ponent of variability in the high-state that may either represent a clear view of the contin- 
uum variability, or else may provide evidence for dense clumps of gas (A^h > 10^^ cm^^) 
in the absorber system [15]. Further observations using Suzaku are key to breaking this 
ambiguity, providing the only method of tracking AGN behavior simultaneously over 
0.5 - 50 keV and the only direct way to probe the importance of Compton-thick gas in 
the absorber complex. 

COMPTON-THICK GAS IN THE LINE-OF-SIGHT 

Suzaku observations of IH 0419-577 and PDS 456 have been key in reshaping our 
view of AGN. For both sources, a marked 'hard excess' of counts was detected in the 
PIN data relative to the predicted flux based on model fits below lOkeV [16, 17]. The 
high PIN-band flux (Fig. 2) can be modeled with a layer of Compton-thick gas (taken 
as Nh — 1.25 X lO^'^cm^^) that covers ~ 90% of the continuum (after correction for 
scattering losses). What is remarkable is that this Compton-thick gas exists in the 
line-of-sight to the nucleus. The type 1 classification of IH 0419-577 implies that the 
Compton-thick gas lies within the BLR, or that the BLR is observed through holes in the 
absorbing medium. Both possibilities challenge current models for AGN and necessitate 
a reconsideration of the true intrinsic luminosity distribution of the population. 

These observations prompted a follow-up study. The Swift BAT survey provides a 
survey of the X-ray sky that is unbiased up to column densities A^h ~ lO^'^cm"^. Tatum 
et al. (in prep) selected a sample of type 1 AGN from the 58-month BAT catalog. To 
understand the sample, we required simultaneous medium and hard-band X-ray data, 
and so the BAT list was cross-correlated with the Suzaku archive, yielding a sample 
comprising 65 observations of 50 objects. The energy flux was measured in the 15- 
50keV and 2-lOkeV bands using the simultaneous XIS and PIN data. The results are 
shown in Fig. 3. The weighted mean ratio is shown as a black solid line; surprisingly, 
1H0419-577 is consistent with the sample mean, having a hardness ratio ~ 1 .4. The red 
line represents the sum of a powerlaw (F = 2. 1) plus reflection from a thin disk of neutral 
material subtending 2k steradians to the illuminating continuum. To characterize the 
reflector we used the pexrav model in XSPEC, assuming no cut-off of the continuum, 
an inclination angle of 60° and solar abundances. The distribution of measured ratios are 
generally much harder than that model. 

The solid orange line represents the hardness ratio for the same reflector, this time 
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FIGURE 2. Suzaku data for 1H0419-577 from 2007(red) and 2010 (green). The solid dark blue line 
shows the total model fit for each epoch. Pale blue lines show the model components, representing three 
different sight-lines through a complex absorber. Aside from the normalization, the spectral model is the 
same for 2007 and 2010 data, and so only component lines from 2010 are shown, for clarity. Note a slight 
offset of the model lines in the PIN band, owing to an 18% re-normalization factor that must be applied 
using the most recent calibration. 



assuming the illuminating continuum is hidden from view (pure reflection). While a pure 
reflection spectrum might be expected for type 2 AGN, such an explanation is at odds 
with a sample of type 1 AGN given the current model for these systems. Markers on the 
right side of the plot represent the ratios expected for the case where a neutral column 
of Compton-thick gas partially covers the continuum source. Several covering fractions 
are denoted: 98% (green), 90% (dark blue), 70% (light blue) and 50% (magenta). In 
the context of this simple partial-covering model, a remarkable 70% of the sample are 
found to have significant covering by Compton-thick gas in the line-of-sight. While 
some Compton-thick AGN have been identified in previous BAT-selected samples of 
AGN [18], our approach, using only the simultaneous constraints provided by the Suzaku 
XIS and PIN, has shown the enormous importance of partially-covering Compton-thick 
gas in the local population of type 1 AGN. 

It is clear that studies below 10 keV have greatly underestimated the intrinsic lumi- 
nosity of a large fraction of the local AGN population (Tatum et al in prep.) and that 
a reconsideration of the structure of the inner regions of AGN is unavoidable. Results 
such as this illustrate the enormous value the broad-band coverage of Suzaku offers 
to the astronomical community. 

BROAD COMPONENTS OF FE Ka EMISSION 

Compelling evidence for Compton-thick gas in the line-of-sight, combined with the cer- 
tainty that the X-ray absorber is outflowing, has motivated detailed calculation of the 
signatures from a Compton-thick flow [19, 3]. Consideration of high wind densities pre- 
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FIGURE 3. Suzaku hardness ratio plotted against BAT flux, using the simultaneous energy density 
fluxes measured from Suzaku PIN/XIS. The sample comprises type 1 AGN selected from the 58-month 
BAT catalog cross-correlated with the Suzaku archive (Tatum et al in prep). The black line represents the 
weighted mean hardness for the sample; the red line represents the hardness ratio expected from a power- 
law reflected from a standard thin disk subtending 2n steradians on the sky; the orange line is the ratio 
expected from pure, neutral reflection. The green, dark blue, light blue and magenta markers represent 
the hardness ratios expected from a partially-covered power-law, with neutral gas and covering fractions 
98%, 90%, 70% and 50%. The continuum power-law is taken to be F = 2.1, for all models. 



diets a strong Fe K emission line should be observed from the wind. This line is broad- 
ened by a combination of electron scattering and velocity effects. Calculations assume 
the wind is launched beyond ~ 30 rg and so relativistic blurring does not contribute to 
the predicted line width, indeed the predicted line widths are not extreme, having typical 
values FWHM~ 1 keV. Case studies of Mrk 766 [19] and PG121 1 + 143 [3] showed that 
a Compton-thick wind can describe the Fe xxv and Fe XXVI absorption lines and the 
broad Fe emission component with a single flow. 

To follow up on those results, Tatum et al (in prep.) have taken an exploratory sample 
of local AGN [see 20] that possess moderately-broad components of Fe K emission. 
Suzaku XIS data were fit over 2-10 keV using a powerlaw continuum modified by a disk 
wind, plus Galactic absorption. Fig. 4 shows some examples of the results obtained. 
The solid blue line represents the continuum power-law convolved with absorption and 
scattering effects from the wind (convolved with the effect of Galactic absorption). The 
green line shows the narrow component of FeKa emission that was required in all 
sources, likely representing reprocessing in distant material, not necessarily associated 
with the wind. The wind model provided a satisfactory fit to all of the broad lines 
observed, indicating that winds may explain both the dominant absorption and emission 
properties of Seyfert type AGN (Tatum et al in prep.). 
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FIGURE 4. Suzaku XIS data for a sample of unobscured AGN. Left to right: Ark 120, Fairall-9, MCG- 
02-14-009 and Swift2127.4 + 5654. The solid blue line shows the Compton-thick wind component of the 
model, the green line is a narrow emission component (see text for details). 



TIMING ANALYSIS 

Measurement of reverberation signals is a powerful technique for constraining the 
size-scale of astronomical systems. Reverberation signatures are measured by cross- 
correlating two time series and searching for time delays between them. To perform this 
technique within the X-ray band, one must understand the contributions of direct and 
reprocessed signatures in each of the bands considered. 

Miller et al. [6, 21] have developed a maximum-likelihood approach to investigate 
time delay signatures in Fourier space. This method finds the power-spectra and cross- 
power-spectra that best fit the time-domain data, with full consideration of data sampling 
and noise, and using estimated measurement uncertainties that account for the variance 
of the source from observation to observation. 

A key part of the timing analysis consists of construction of a "lag spectrum" that 
shows the time lags between two time series as a function of the frequency of the source 
variation. Analysis of the lag spectrum yields information about the geometry and size 
scale of the reprocessor. Consider a thin spherical shell of material surrounding an active 
nucleus, where the directly-viewed continuum is emitted from a central point source 
and measured in one time series, while light scattered from the shell is measured in 
a second time series. If the shell has holes "partial-covering" in the sight-line to the 
observer, and also on the opposite side of the continuum source from the observer, then 
the gas appears somewhat like a ring seen end-on to the observer. One can construct the 
lag spectrum for such a reprocessor. At high frequencies, where the variability mode is 
significantly shorter than the light travel time across the shell, scattering from different 
parts of the shell add incoherently over a wide range of phases, so the observed net time 
delay tends to zero. One also observes some oscillations that depend on the details of the 
shell structure. At low frequencies, where the source variability period is much longer 
than the light travel time across the shell, the scattered light adds mostly coherently 
and a mean time delay is observed for the shell. At a frequency where the time period 
equals the light travel time across the shell, a sharp transition is observed between these 
two regimes and the measured time delay may become negative for some ranges of 




FIGURE 5. Lag spectra between hard and soft bands in X-ray observations. Left to right: Suzaku 
data from NGC4051 [6]; XMM-Newton data from MCG-6-30-15 (Miller et al., in prep; also see Em- 
manoulopoulos et al. 22), XMM-Newton data from 1H0707^95 [21]. Lags are positive when the hard 
band lags the soft band. 



frequency. While all the scattered light signals actually are delayed, the negative lag 
appears in Fourier space if the phases wrap around by more than Itu radians. Thus, the 
detailed oscillatory structure in the lag spectrum tells us about the structure, geometry 
and size- scale of the reprocessor. Note that if the shell were symmetric and full-covering, 
the time delays oscillate but stay positive. 

Fig. 5 shows some examples of lag spectra having clear transitions from positive lags 
(hard band lagging soft) at low frequencies to zero or negative lags at high frequencies. 
The observed transitions are broader than expected for a thin-shell reprocessor, indicat- 
ing that the reverberating material has significant depth [6]. The observed time delays 
have values of hundreds of seconds at low variation frequencies but fall to values close 
to zero at high frequencies. A typical transition frequency is V ~ 5 x 10~^Hz, indicating 
the reprocessor to be a few light hours from the central source, equivalent to a fews tens 
to a few hundreds of gravitational radii, depending on the black hole mass. 

The interpretation of large-scale X-ray reverberation is supported by a substantial 
body of spectral evidence [23]. Zones of partially-covering, absorbing circumnuclear 
material shape the X-ray timing and spectral properties of AGN. Both spectral and 
timing analyses indicate that the global covering factor of the reprocessor must be 
high (~ 50%) in order to achieve a sufficient amount of scattered light to explain the 
observations. Further work is needed to create models that reproduce both the timing 
behaviour and the observed high resolution spectra. 
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